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Quantum superconductor-insulator transition: Implications of BKT-critical behavior 
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We explore the implications of Berezinskii-Kosterlitz-Thouless (BKT) critical behavior on the 
two dimensional (2D) quantum superconductor-insulator (QSI) transition driven by the tuning 
parameter x. Concentrating on the sheet resistance R (x, T) BKT behavior implies: an explicit 
quantum scaling function for R(x,T) along the superconducting branch ending at the nonuniversal 
critical value R c = R(x c ); a BKT-transition line T c (x) oc (x — x^)" v where z is the dynamic and 
V the exponent of the zero temperature correlation length; independent estimates of zv, z and V 
from the x dependence of the nonuniversal parameters entering the BKT expression for the sheet 
resistance. To illustrate the potential and the implications of this scenario we analyze data of 
Bollinger et al. taken on gate voltage tuned epitaxial films of La2-xSr x Cu04 that are one unit 
cell thick. The resulting estimates z ~ 2.35 and v ~ 0.63 point to a 2D-QSI critical point where 
hyperscaling, the proportionality between d/X 2 (0) and T c , and the correspondence between quantum 
phase transitions in D and classical ones in (D+z) dimensions are violated and disorder is relevant. 

PACS numbers: 74.40. Kb, 74.62.-c, 74.78.-w 



I. INTRODUCTION 

A variety of different materials undergo a quantum 
superconductor- insulator (QSI) transition in the limit 
of two dimensions (2D) and zero temperature by vari- 
ation of a tuning parameter including film thickness, 
disorder, applied magnetic field, and gate voltage^— 
A widespread observable to study this behavior is the 
temperature dependence of the sheet resistance R (x, T) 
taken at various values of the tuning parameter x. The 
curves of R(T) at different x resemble the flow to a 
nearly temperature independent separatrix between su- 
perconducting and insulating phase with sheet resistance 
R c = R (x = x c , T ~ 0). This behavior implies a cross- 
ing point of the isotherms R (x) at different tempera- 
tures at x c which is a characteristic feature of a quantum 
phase transition and in the present case of a QSI tran- 
sition. Traditionally the interpretation of experimental 
data taken close to the 2D-QSI transition is based on 
the quantum scaling relation R (x, T) = R c (x) G (y) with 
y = c(x — x c ) /T x / Zu z is the dynamic, v the critical 
exponent of the zero temperature correlation length and 
c a nonuniversal coefficient of proportionality^— Given 
sheet resistance data fits to this scaling form yield esti- 
mates for the critical value of the tuning parameter x c 
and the exponent product zv, properties which are in- 
sufficient to distinguish different models, to fix the uni- 
versality class to which the QSI transition belongs, or to 
clarify the relevance of disorder . 

The nature of the 2D-QSI transition has been in- 
tensely debated^— The scenarios can be grouped into 
two classes, fermionic and bosonic. In the fermionic case 
the reduction of T c and the magnitude of the order pa- 
rameter is attributed to a combination of reduced density 
of states, enhanced Coulomb interaction and depairing 
due to an increase of the inelastic electron-electron scat- 
tering ratej*^ The bosonic approach assumes that the 
fermionic degrees of freedom can be integrated out, the 
mean square of the order parameter does not vanish at 



T c , phase fluctuations dominate and the reduction of T c 
is attributable to quantum fluctuations and in disordered 
systems to randomness in additionj&iii This scenario is 
closely related to the suppression of ferroelectricity^ e.g. 
in SrTi0 3 3 

Here we adopt the bosonic scenario and concentrate 
on systems which undergo at finite temperature a nor- 
mal state to superconductor transition with Bcrczinskii- 
Kosterlitz-Thouless (BKT) critical behavior ; 14 ' 15 origi- 
nally derived for the 2D xy-modcl with a two component 
order parameter and short range interactions. The oc- 
currence of BKT criticality in 2D superconductors also 
implies that the mean square of the order parameter does 
not vanish at T c and with that there are, in analogy to 
4 Hc, condensed pairs (bosons) below and uncondensed 
ones above T c . In 4 He and superconductors the order 
parameter is a complex scalar corresponding to the com- 
ponents in the xy-model. Supposing that in supercon- 
ductors the interaction of Cooper pairs is short ranged 
and their effective charge is sufficiently small the criti- 
cal properties at finite temperature are then those of the 
3D-xy (bulk) and 2D-xy (thin films) models ) 16 ' 17 rem- 
iniscent to the lamda transition in bulk 4 Hei£ and the 
BKT-transition in thin 4 He filmsi 1 ^— In this context 
it is important to recognize that the existence of the 
BKT-transition (vortex-antivortex dissociation instabil- 
ity) in 4 He films is intimately connected with the fact 
that the interaction energy between vortex pairs depends 
logarithmic on the separation between them. As shown 
by Pearl ^ vortex pairs in thin superconducting films 
(charged supcrfluid) have a logarithmic interaction en- 
ergy out to the characteristic length A2D = A 2 /d, beyond 
which the interaction energy falls off as 1/R. Here A is 
the magnetic penetration depth and d the film thickness. 
As A2D increases by approaching T c the diamagnetism 
of the superconductor becomes less important and the 
vortices in a clean and thin superconducting film become 
progressively like those in 4 He films. 24 ' 25 

The occurrence of a 2D-QSI transition implies 
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a line T c (x) of BKT transition temperatures end- 
ing at the quantum critical point at x = x c 
where T c (x = x c ) — 0. It separates the super- 
conducting from the insulating ground state. The 
BKT transition is rather special because the cor- 
relation length diverges above T c as £(x,T) = 

£o (x) exp (x) /2T C 1/2 (x)) (T/T c (x) - 1)~ 1/2 ) and 

the sheet resistance tends to zero according to R (x, T) = 

R (x) exp (- (b R (x) /T c 1/2 (x)) (T/T c (x) - 1)- 1/2 ) 

Approaching the 2D-QSI transition quantum phase fluc- 
tuations renormalize Rq(x), bn{x), and T c (x). Indeed 
the BKT-transition line approaches the 2D-QSI tran- 
sition as T c = (c (x — x c ) /yc) zu because the quantum 
scaling form G (y) exhibits at the universal value y c of the 
scaling argument a finite temperature singularity] 6 ' 8 ' 16 
Noting that the amplitude of the BKT correlation 
length £o (x) should match the divergence of the quan- 
tum counterpart, £ (T = 0) oc (x — x c )~ oc T~ 1/z the 
exponents z and V should emerge from the amplitude 
Ro (x) in terms of R (x) - R c oc ^ 2 (x) cx (T = 0) cx 

{x — x c ) 21 ' oc T 2 ^ z Accordingly, the nonuniversal func- 
tions T c (x) and Rq (x) entering the BKT expression for 
the sheet resistance exhibit close to the QSI transition 
quantum critical properties disclosing the quantum 
critical exponents zV, z and V. The exponents z and 
V are characteristic properties of the universality class 
to which the QSI transition belongs. In addition, given 
their values the relevance of disorder and the equivalence 
between quantum phase transitions in systems with 
D spatial dimensions and the ones of classical phase 
transitions in (D + z) dimensions can be checked£&i£ 
The occurrence of a BKT transition line also implies: 
a nonuniversal critical sheet resistance R c — Rq (x c ) 
because Rq (x) is nonuniversal; an explicit form of the 
superconductor branch of the quantum scaling function 
G(c (x-xJ/T 1 /™). 

Even though BKT critical behavior is not affected by 
short-range correlated and uncorrelated disorde r 27 ! 28 the 
observation of this behavior requires sufficiently homoge- 
neous films and a tuning parameter which does not affect 
the disorder. Noting that sample inhomogeneity and vor- 
tex pinning are relevant in thickness and perpendicular 
magnetic field tuned transitions, electrostatic tuning of 
the 2D-QSI transition using the electric field effect ap- 
pears to be more promising] 26 ! 29 " — Indeed electrostatic 
tuning is not expected to alter physical or chemical dis- 
order, but changes the mobile carrier density. 

In Sec. II we sketch the theoretical background. To il- 
lustrate the potential and the implications of finite tem- 
perature BKT criticality on the 2D-QSI transition we 
analyze in Sec. II the sheet resistance data of Bollinger 
et al£2- taken on gate voltage tuned epitaxial films of 
La2- 2; Sr :r Cu04 that are one unit cell thick. Commenting 
on the difficulties in observing the BKT features in the 
magnetic penetration depth we consider the data of Bert 
et al^ taken on the superconducting LaA10a/SrTi03 



interface. 



II. THEORETICAL BACKGROUND 

Continuous quantum-phase transitions (QPT) are 
transitions at zero temperature in which the ground 
state of a system is changed by varying a parameter 
of the Hamilt onian £ ' 1 6 ' 1 7 The quantum superconductor- 
insulator transitions (QSI) in two-dimensional (2D) sys- 
tems tuned by disorder, film thickness, magnetic field or 
with the electrostatic field effect are believed to be such 
transitions.— ~— ' 16 ' 17 ' 26 ' 30 Traditionally the interpretation 
of experimental data taken close to the 2D-QSI transi- 
tion is based on the quantum scaling relation^ - — 



R{x,T) 
R, 



G{y),y 



c \x 



J'l/zV 



(1) 



where R is the resistance per square, R c the limiting 
(x — > x c and T — > 0) resistance, x denotes the tuning 
parameter and c a nonuniversal coefficient of proportion- 
ality. G (y) is a universal scaling function of its argument 
such that G (y = 0) = 1. In addition, z is the dynamic 
and v the critical exponent of the correlation length, sup- 
posed to diverge as 



Z(T = 0)=£ o \x-x c \- 



(2) 



The critical sheet resistance R c separating the supercon- 
ducting and insulating ground states is determined from 
the isothermal sheet resistance at the crossing point in 
R (T) vs. tuning parameter x at x c . The existence of such 
a crossing point, remaining temperature independent in 
the zero temperature limit, is the signature of a QPT. 
The data for R (x, T) plotted vs. \x - x c \ /T x / Vz should 
then collapse onto two branches joining at R c . The lower 
branch stems from the superconducting {x — x c > 0) and 
the upper one from the insulating phase (x — x c < 0). 
This scaling form follows by noting that the divergence 
of the zero temperature correlation length, £ (T = 0) = 
£o \5\ is at finite temperature limited by the length 
Lt oc T~ x l z ^ Thus G (y) is a finite-size scaling function 
because y oc [L T /£ (T = 0)] 1/T7 oc \x-x c \ jT x l Vz . Sup- 
posing that there is a line of finite temperature phase 
transitions T c (x) ending at the quantum critical point 
T c (x = x c ) = the quantum scaling form (TTJ) exhibits 
at the universal value y c of the scaling argument a finite 
temperature singularity] 6 ' 8 ' 16 The phase transition line is 
then fixed by 



c\x - x c \ 



lie 



(3) 



Otherwise one expects that sufficiently homogeneous 2D 
superconductors exhibit at the superconductor to normal 
state transition BKT critical behavior i 14 ' 15 ' 26 Note that 
there is the Harris criterion^ 7 - stating that short-range 
correlated and uncorrelated disorder is irrelevant at the 
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unperturbed critical point, provided that v > 2/D, where 
D is the dimensionality of the system and v the criti- 
cal exponent of the finite-temperature correlation length. 
With D = 2 and v = oo , appropriate for the BKT 
transition ; 14 ' 15 this disorder should be irrelevant. Given 
a BKT superconductor to normal state transition the 
sheet resistance scales for T >T C (x) > as^ir— 



R(x,T) 
Ro (a?) 



exp 



bn (x) 



tV 2 {x) (T/T c (x)-1) 



1/2 



(4) 



allowing to probe the characteristic BKT correlation 
length 14 ) 15 ' 25 . 26 



gx,T) 

Zo(x) 



exp 



bR (x) 



2T C 1/2 (x) (T/T c (x) - 1) 



1/2 



(5) 



While Rq(x), bR(x), and T c {x) depend on the tun- 
ing parameter and are subject to quantum fluctuations 
the characteristic BKT form of the correlation length 
and with that of the sheet resistance applies for any 
T >T C (x) > 0. Through standard arguments (see, e.g., 
Refi^) quantum mechanics does not modify universal fi- 
nite temperature properties. bR (x) is given b y 15 ' 26 



bR (x) = b R T t 



-1/2 



(0) 



with bR = Ait jb. The parameter b is expected to remain 
constant in the low T c regime.— It is related to the vor- 
tex core energy E c in terms of b = f (E c /kBT c )r^ and 
controls below the Nelson-Kosterlitz jump the temper- 
ature dependence of the magnetic penetration depth in 
terms of (A (T c ) /A {T)) 2 = l + (b/4) (T/T c - 1) 1/2 .^ The 
amplitude of the BKT correlation length £ (x) is propor- 
tional to the vortex core radius^ known to increase with 
reduced T c 2&£- Indeed, the zero temperature correlation 
length i (T = 0) diverges as £ (T = 0) oc (x - x c )~ v (Eq. 
© and combined with T c oc {x - xc) zv (Eq. ©) we ob- 
tain £ (T = 0) oc (x - x c )~ v oc Tc 1/z . Noting that R (x) 
approaches R c from above the scaling relation 

R (x) - R c oc r 2 (T = 0) oc (x - x c ) W oc T 2 J\ (7) 

should apply— making the determination of the expo- 
nents V and z possible. Consistency requires that the 
resulting zv agrees with the estimate obtained from the 
critical BKT line T c (x) oc (x — x c ) zu . Other implications 
concern the universality class of the 2D-QSI transition 
and the relevance of disorder. Given estimates for V and 
z the equivalence between quantum phase transitions in 
clean systems with D spatial dimensions and the ones of 
classical phase transitions in {D + z) dimensions can be 
checked. The fate of a clean critical point under the influ- 
ence of disorder is controlled by the Harris criterion; 27 ' 28 
If the zero temperature correlation length critical expo- 
nent fulfills the inequality 17 > 2/D the disorder does not 
affect the critical behavior. If the Harris criterion is vio- 
lated, v < 2/D, the generic result is a new critical point 



with conventional power law scaling but new exponents 
which fulfill the Harris criterion. Another option is that 
the disorder destroys the QSI transition. 

Finally, given a BKT transition line expression (U) the 
sheet resistance given by Eq. ([5]) transforms with Eq. (]6]) 
and the scaling variable y = c \x — x c \ /T x l zv to 



In 



R(x,T) 
Ro (x) 




-1/2 



(8) 



valid foe y < y c . Close to quantum criticality where 
R (x) ~ R c reduces to 



In 



R{x,T) 
R c 




-1/2 



- 1 



= WGfo)). 



(9) 



These relations are explicit forms of the quantum scal- 
ing function G (y) applicable to the superconductor 
branch. They reveal that the critical sheet resistance 
Rq (x —> x c ) = R c is the endpoint of a nonuniversal func- 
tion and accordingly nonuniversal. Another implication 
concerns the universality class of the 2D-QSI transition. 
Supposing that the equivalence between quantum phase 
transitions in clean systems with D spatial dimensions 
and the ones of classical phase transitions in (D + z) di- 
mensions applies, the 2D-QSI transition at the endpoint 
of a BKT line T c (x) should belong to the finite tempera- 
ture (2 + z) — xy universality class, xy denotes an order 
parameter with two components, including the complex 
scalar, ^ = |$| exp (i(p), of a superconductor j 5 -* 1 ^ 

The BKT - theory of thermally-excited vortex- 
antivortex pairs also predicts a super-to-normal state 
phase transition marked by the Nelson-Kosterlitz jump, 
a discontinuous drop in supcrfluid density frorn^ 5 -, 



A 2 Tc 



32n 2 k B T c 



1.02T C , 



(10) 



to zero. The numerical relationship applies for c?/A 2 (T~) 
in cm -1 and T c in K. d denotes the thickness of the 2D 
system, A the in-plane magnetic penetration depth, and 
<l>o = hc/2e. In addition there is the prediction that 
d/X 2 (T = 0), a measure of the phase stiffness, scales near 
the endpoint of the BKT transition line a o 6 ' 8 ' 16 ' 17 



A 2 (T = 0) 



167T 3 fc B T c 



Q 2 1.6Q 2 T C , 



(11) 



provided that D + z is below the upper critical dimension 
D u where hyperscaling holds. Since D u = 4 in the D- 
xy-model the validity of Eq. (fTTj) is in D = 2 restricted 
to z < 2. Relation (fTTj) is the quantum counterpart of 
the Nelson-Kosterlitz relation ©. Q2 is a dimensionless 
critical amplitude with the lower boun d 16 ' 17 



Qi > 2/tt, 



112) 
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dictated by the characteristic temperature dependence 
of d/X 2 (T) below the Nelson-Kosterlitz jump (Eq. 
(pjD.iL 1 ! Combining Eqs. CLO} and (JTTJ we obtain 

where p s is the supcrfluid density. The supcrfluid transi- 
tion temperature T c as a function of the supcrfluid den- 
sity p s (T = 0) has been measured in 4 He films for tran- 
sition temperatures ranging from 0.3 to 1 K by Crowell 
et al£- They studied 4 He films adsorbed in two porous 
glasses, aerogel and Vycor, using high-precision torsional 
oscillator and dc calorimetry techniques. The inves- 
tigation focused on the onset of superfluidity at low 
temperatures as the 4 He coverage is increased. Their 
data yields p s (T = 0) ~ 15.3T C with p s in ^moles/m 2 
and T c in K. Combined with the BKT transition line 
p s (T c = 0) = 8.73T C we obtain Q 2 ~ 1.12. 



Comparing the data with the respective lines we observe 
that the BKT regime is attained and that hat the BKT 
T c is almost an order of magnitude lower than the mean- 
field counterpart. This uncovers a BKT transition from 
uncondensed to condensed Cooper pairs driven by strong 
phase fluctuations. On the other hand it is important to 
recognize that agreement with BKT-criticality is estab- 
lished in a temperature window only. Its upper bound 
reflects the crossover from BKT- to AL-behavior while 
the lower bound stems from the rounded BKT-transition. 
Precursors of this phenomenon are clearly visible in Fig. 
[T] around 7 K. Here the correlation length is prevented 
to grow beyond a limiting length L, i.e. the linear ex- 
tent of the homogeneous domains. As a result, a finite 
size effect and with that a rounded transition occurs^ 
Because the BKT correlation length does not exhibit the 
usual power law divergence of the correlation length as 
T c is approached, it is particularly susceptible to such 
finite-size effect. 



III. COMPARISON WITH EXPERIMENT 

To illustrate the potential and the implications 
of the outlined BKT scenario we analyze next the 
data of Bollinger et ai— taken on epitaxial films of 
La2_ a; Sr a ;Cu04 that are one unit cell thick. Very large 
electric fields and the associated changes in surface car- 
rier density enabled shifts in the midpoint transition tem- 
perature T c by up to 30 K. Hundreds of resistance ver- 
sus temperature and carrier density curves were recorded 
and shown to collapse onto a single function, as the 
quantum scaling form (Eq. (p}) for a 2D-QSI transi- 
tion predicts. The observed critical resistance is close to 
the quantum resistance for pairs, Rq = h/4e 2 ~ 6.45 
kS7. Our starting point is the temperature dependence 
of the sheet resistance taken at various gate voltages V g 
where a BKT transition is expected to occur. As an ex- 
ample we depicted in Fig. [T]a, the data for V g = —1.5 
V. The observation of BKT-behavior requires that the 
data extend considerably below the mean-field transi- 
tion temperature T c0 . We estimated it with the aid 
of the Aslamosov-Larkin (AL) expression^ for the con- 
ductivity a{T,V g ) = a n (V g ) + a / (T/T c0 - 1), with 
(T = ne 2 /8h ~ 1.52 x 10~ 5 fi -1 , where Gaussian fluctu- 
ations are taken into account and T c q is the mean-field 
transition temperature. The resulting temperature de- 
pendence is included in Fig. QJi. It clearly reveals that 
the data extend considerably below T c q ~ 12.5 K. To 
establish and characterize BKT behavior below T c q we 
invoke Eq. (0} and 

As indicated in Fig. QJ> this relation is used to fix &^ (T c ) 
and T c (x) while R (x) is estimated by adjusting Eq. (3]) 
with given bn (T c ) and T c (x) to the sheet resistance data. 




T(K) 

FIG. 1: (color online) a) Sheet resistance R(T) for V g — —1.5 
V taken from Bollinger et a/.— The dotted curve is a fit to the 
Aslamosov-Larkin expression yielding the T c o — 12.5 K. The 
solid line is the BKT behavior Eq. (@} with Ro = 18 kfi and 
the T c , b R values derived in Fig. [lb. b) (din (R) /dT)~ 2/3 vs. 
T derived from the data shown in Fig. [T^. The dotted line 
indicates the AL-behavior and the solid one Eq. (| 13p with 
T c = 1.8 K and (2/b R (T c )) 2/3 = 0.49 (K~ 1/3 ) 
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To explore the quantum critical behavior latent in 
T c (V g ), R (V g ) and 6r (T c ) we performed this analysis 
of the temperature dependence of the sheet resistance 
for additional gate voltages. To fix the critical gate volt- 
age V gc we used the empirical gate voltage dependence 
of the number of mobile holes x per one formula unit of 
Bollinger et alM yielding x (V g ) = x c + 0.012(Vg C - V g ) 
down to V g = -2V with x c = 0.0605 and V gc ~ -0.7 
V. The results T c (V g ), R (V g ) are shown in Fig. [2ji and 
Fig. [5b while bn (T c ) is depicted in Fig. [3] As can be 
seen in Fig. [2Jt the BKT transition line differs substan- 
tially from the so called superconducting dome behavior 
observed in bulk cuprate superconductors. It is approxi- 
mately given by T c (x) /T C max = 1 — 82.6 (x — x m ) where 
r c max(ira =0.16) is the maximum T c £^ Close to the 
QSI transition where even bulk cuprate superconductors 
become essentially 2Di£ it reduces to T c (x) /T c max ~ 
10.3 (x — x c ) with x c ~ 0.049 and suggests zV ~ 1. In 
any case it differs substantially from the BKT line shown 
in Fig. [5^ yielding the estimate 

zv ~ 1.46, (15) 

in agreement with the value zv ~ 1.5 derived by Bollinger 
et al— using the quantum scaling approach. In contrast 
to this, from the nonuniversal parameters Rq (V g ) and 
Ro (T c ) shown in Fig. [2Jd we derive with Eq. (O in 
addition 

z7~ 0.63,2 ~ 2.35, zv= 1.48. (16) 

As these exponents satisfy the inequality (2 + z) v > 2 we 
use the correct scaling argument, because 6 = (fj, — /i c ) oc 
(x — x c ) oc (Vgc — Vg) where /i denotes the schemical 
potential^ The agreement between these zv values con- 
firms the applicability of the scaling relation ([7|). Simi- 
larly, the nonuniversal parameter bj^ (x) exhibits accord- 
ing to Fig. [3] the expected T c dependence (Eq. ((6|)). 
Noting that D + z ~ 2 + z ~ 4.35 exceeds the upper 
critical dimension D u = 4 the critical exponent of the 
zero temperature correlation length v should adopt its 
mean-field value v = 1/2. However the fate of this clean 
critical point under the influence of disorder is controlled 
by the Harris criterion i 27 i 28 If the inequality 77 > 2/D is 
fulfilled, the disorder does not affect the critical behavior. 
If the Harris criterion is violated ( 77 < 2/D), the generic 
result is a new critical point with conventional power law 
scaling but new exponents which fulfill 77 < 2/D. Since 
77 = 1/2 violates this inequality disorder is relevant and 
drives the system from the mean-field to an other criti- 
cal point with different critical exponents as our estimate 
77 ~ 0.63, fulfilling 77 < 2/D = 1, uncovers. The resulting 
2D-QSI transition with v ~ 0.63 and z ~ 2.35 violates 
then the equivalence between quantum phase transitions 
in systems with D spatial dimensions and the ones of 
classical phase transitions in (D + z) dimensions. In ad- 
dition the proportionality between d/A 2 (0) and T c (Eq. 
(fTTT) ') valid below the upper critical dimension D u = 4 
does no longer hold because (D + z) ~ 4.35 is above 



D u = 4. In fact magnetic penetration depth measure- 
ments taken on underdoped high quality YBa2Cu3C , 6+ :c 

single crystals revealed T c oc (e?/A 2 (0)) ' .— Another 
option is that the disorder destroys the QSI transition. 
Given the evidence for a rounded BKT transition (see 
Fig. [1]) and the missing data close to the QSI transition 
(see Fig. @| further studies are required to elucidate this 
option. 




o-l , 1 , 1 , 1 , 

-2.5 -2.0 -1.5 -1.0 -0.5 

V v > 

FIG. 2: (color online) a) Phase transition line T c (V g ) de- 
rived from the sheet resistance data of Bollinger et a/.— us- 
ing Eq. (|14[) . The solid line is a least square fit yielding 
T c = 2.41 (Vgo - VgY 17 (K) and zu = 1.46 with V gc = -0.7 V. 
b) Ro (V g ) derived from the sheet resistance data of Bollinger 
et alM using Eq. © with given T c (V g ) and b R (T c ); Ro (T c ) 
derived from Rq (V g ) using T c (V g ). The solid line is Ro (V g ) = 
Rc + (V gc - V S ) 2TJ (kfi) with R c = 6.45 kfl, V gc = -0.7 V and 
77 = 0.63. The dashed line is R (T c ) = R c + 7.52T c 2/z (kft) 
with z — 2.35. 

To complete the analysis of the data of Bollinger et 
alM W e depicted in Fig. H the plot R(V g ,T) / R (V g ) 
vs.(T c (V g ) /T) 2 ^ 3 corresponding to the quantum scal- 
ing function G (y) (Eq. ([8])) in terms of y/y c = 
(T c (V g ) /T) 2 ^ 3 . Apparently, the data does not fall com- 
pletely on the BKT curve indicated by the dashed line. 
It corresponds to the superconductor branch of the quan- 
tum scaling function. Instead we observe a flow to and 
away from the universal characteristics. As T c /T de- 
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0-1 , 1 , 1 , 1 , 

12 3 4 

T (K) 



FIG. 3: (color online) Estimates for bn (T c ) derived from the 
sheet resistance data from Bollinger et aL— using Eq. (|13|l . 
The solid line is b R (T c ) = b R Tl' 2 (K 1/2 ) with b R = 4.61 in 
agreement with Eq. ([6]). 

creases for fixed T c the crossover to AL-behavior sets in, 
while the rounding of the transition leads with increas- 
ing T c /T to a flow away from criticality. The important 
lesson then is that the quality of the data collapse on a 
single curve heavily depends on the temperature range 
of the data entering the plot. Another striking feature is 
the extended scaling regime. Within the BKT scenario 
it simply follows from the fact that the scaling form (|4]) 
applies along the BKT transition line irrespective of the 
distance from the QSI transition. In the quantum scal- 
ing approach this property remains hidden and merely 
suggests an extended quantum critical regime. The ex- 
cellent quality of the piecewise data collapse also reveals 
that the observance of the substantial variation of Ro (V g ) 
(sec Fig. [5Jd) is essential, while in the quantum scaling 
approach it is fixed by the critical sheet resistance. To 
demonstrate these features even more compelling we de- 
picted in Fig. [Sh, R/R c vs.(V gc - V g ) /T 2 / 3 . Apparently 
the data do not fall even piecewise on a single curve. 
After all this is not surprising because the quantum scal- 
ing form ([T]) holds close to the critical sheet resistance 
only and Ro (V g ) varies substantially in the gate voltage 
regime considered here (see Fig. [2t>). As shown in Fig. 
[5)d and Fig. [6] this behavior allows to estimate Ro (V g ) 
from R/R c vs.(V gc — V g ) /T 2 / 3 by rescaling R c in terms 
of b (V g ) R c with b (V g = —IV) = 1 to achieve piecewise a 
collapse of the data. The resulting b (V g ) is shown in Fig. 
[6] and agrees well with Ro(V g )/Ro(V g = —IV) derived 
from the BKT behavior of the sheet resistance (see Fig. 

A BKT line with a QSI transition at its endpoint was 
also explored rather detailed at the interface between 
the insulating oxides LaA103 and SrTiC>3 exhibiting a 
superconducting 2D electron system that can be mod- 
ulated by a gate voltagei 26 ' 32 ' 39 ' 40 BKT behavior and 
with that a 2D electron system was established as fol- 
lows: The the current-voltage characteristics 3 -^ revealed 
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FIG. 4: (color online) R(V g ,T) /R (V g ) vs.(T c (V g ) /T) 2/A = 
y/y c corresponding to the BKT quantum scaling function 
G(y) (Eq. ©). R(V g ,T) is taken from Bollinger et alM 
in the temperature range from 4.4 K to 100 K. The respec- 
tive values for Ro (x), T c (V g ) and b r are shown in Figs. [2] 
and [3] The dashed line is the critical BKT behavior given by 
Eq. © with b R = b R /Tc /2 = 4.55. 



at the BKT transition temperature T c the characteris- 
tic BKT form V oc I a with a = 3.— Consistency with 
the characteristic temperature dependence of the sheet 
resistance (Eq. $$) was established! 26 ' 39 ' 40 It was also 
shown that the effective thickness of the superconduct- 
ing 2D system can be extracted from the magnetic field 
dependence of the conductivity at T c 4i The gate volt- 
age tuned BKT phase transition line, T c (V g ), derived 
from the temperature dependence of the sheet resistance 
at various gate voltages uncovered with Eq. ((4| consis- 
tency with T c (V g ) = 8.9 x 1(T 3 (V g - V gc ) 2/3 K pointing 
to quantum critical behavior (Eq. ([3])) with zv = 2/3 
and the critical sheet resistance R c = 2.7 kf h 26 ' 40 Fur- 
thermore the estimates 77 ~ 2/3 and z ~ 1, confirming 
zV = 2/3, have been derived from Ro (V^).— This sug- 
gests that the gate voltage tuned QSI transition of the 2D 
electron system at the LaAlOs/ SrTiC>3 interface belongs 
to the (2 + z) = 3-xy universality class where hyperscal- 
ing and with that the proportionality between d/X 2 (0) 
and T c (Eq. (fTTj) ) applies. On the other hand because 
77 ~ 2/3 < 2 /D = 1 disorder is according to Harris theo- 
rem relevant as well . 27 ' 28 

To comment the difficulties in observing the BKT fea- 
tures in the magnetic penetration depth as well as the 
relation between d/X 2 (0) and T c , we reproduced in Fig. 
[7] the data of Bert et alS^ for a gate voltage tuned su- 
perconducting LaAlOa/SrTiOa interface in terms of T c 
vs. d/X 2 (T = 0.0AK). T c is defined here as the tem- 
perature at which the diamagnetic screening drops be- 
low the noise level corresponding to a detectable d/X 2 of 
0.10 - 0.34 cm -1 . A glance at Fig. [7] reveals that this 
is just the regime where the universal quantum behavior 
applies, indicated by the dashed and solid lines, corre- 
sponding to the lower bound (|12[) and the behavior de- 
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FIG. 5: (color online) a) 'R/R c vs.(V gc - V g ) /T 2/3 with 
R c = 6.45 kf2 and V gc = —0.7 V derived from the tempera- 
ture dependence of the sheet resistance of Bollinger et a/.— 
The data cover the range from 4 K to 100 K. b) R/b (V g ) R c 
vs.(V gc - V g )/T 2/Z with b(V g = -lV) = 1 and b(V g ) ad- 
justed to achieve a piecewise collapse of the data. 
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FIG. 6: (color online) b(V g ) and R (V g )/R (V g = 
-IV). The solid line is R {V g )/R {V g = -IV) = 
(i? c + 15.9 (V 9C - V,) 1 - 26 ) /R (V g = -IV) with R c = 6.45 
kfi, V m 
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rived from the 4 He data of Crowell et a/.,— respectively. 
Nevertheless the data reveal the flow to the 2D-QSI tran- 
sition which is attained at much lower T c 's. Otherwise 
the data points resemble the outline of a fly's wingjii re- 
markably similar to the T c vs. 1/A^ b (T = 0) plots of the 
bulk superconductors Y . 8 Cu .2-123, Tl-1212^ and Tl- 
2201—, covering nearly the doping regime of the so called 
superconducting dome extending from the underdoped 
to the overdoped limit . According to the generic plot 
T c /T c (x m ) vs. j(x rn ) /~f(T c ), where 7 = £ a b/£ c is the 
anisotropy and £ a b,c denote the in-plane and c-axis cor- 
relation length, these cuprates become nearly 2D in the 
underdoped limit In any case, Fig. 0shows that in this 
plot the universal QSI behavior is attained at compara- 
tively low T c values only. Accordingly, in the QSI regime 
of interest the Nelson-Kosterlitz jump given by Eq. (fTU)) 
becomes very small and appears to be beyond present 
experimental resolution.— On the contrary in the tem- 
perature dependence of the sheet resistance is the BKT 
critical regime accessible because T c /T c q ~ 0.75^ even 
though small compared to that in the La2- s Sr^Cu04 



FIG. 7: (color online) T c vs. d/\ 2 (T = 0.04K) for a gate volt- 
age tuned superconducting LaAlOs/SrTiOs interface taken 
from Bert et al^ The solid and dashed lines are the pre- 
dicted universal relationship © close to the QSI transition 
with Q2 = 2/-7T, resulting from the lower bound (|12[1 . and 
Q2 = 1.12 derived from the 4 He film data of Crowell et a/.— 

films where T c /T c0 w 0.1. 

IV. SUMMARY AND DISCUSSION 

In sum, we sketched and explored the implications 
of Bcrczinskii-Kosterlitz-Thoulcss (BKT) critical behav- 
ior on the quantum critical properties of a two dimen- 
sional (2D) quantum superconductor-insulator transition 
(QSI) driven by the tuning parameter x. It was shown 
that the finite temperature BKT scenario, implies in 
terms of the characteristic temperature dependence of 
the BKT correlation length an explicit quantum scaling 
function for the sheet resistance R(x,T) along the su- 
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perconducting branch ending at the nonuniversal critical 
value R c = Rq (x c ). This scaling form fixes the BKT- 
transition line T c (x) and provides estimates for the quan- 
tum critical exponent product zv. In addition, indepen- 
dent estimates of z77, z and 77 follow from the x depen- 
dence of the nonuniversal parameters T c (x) and Rq (x) 
entering the characteristic BKT expression for the sheet 
resistance R(x,T). This requires that the BKT criti- 
cal regime where phase fluctuations dominate is attained 
and the finite temperature BKT relation for the sheet 
resistance applies for any T >T C (x) > 0. The last con- 
dition is satisfied because the BKT expression for the 
sheet resistance is simply related to the characteristic 
temperature dependence of the BKT correlation length. 
Quantum fluctuations enter via the nonuniversal parame- 
ters T c (x) and Rq (x) disclosing close to x c the respective 
quantum critical behavior. Even though BKT critical be- 
havior is not affected by short-range correlated and un- 
corrected disorde r 27 ' 28 the observation of this requires 
sufficiently homogeneous films and a tuning parameter 
which does not affect the disorder. Noting that in the 
magnetic field tuned case there is no BKT-line, thick- 
ness and gate voltage tuned 2D-QSI transitions appear 
to be promising candidates. In any case the scenario out- 
lined here requires a line of BKT-transitions T c (x) with 
a quantum critical endpoint T c (x c ) = and sheet resis- 
tance data which attain the BKT critical regime. 

To illustrate the potential and the implications of this 
scenario we analyzed data of Bollinger et al£2- taken on 
gate voltage tuned epitaxial films of I^-^Sr-jCuCU that 
are one unit cell thick. Evidence for dominant phase 
fluctuations and BKT-critical behavior was established 
in terms of the temperature dependence of the sheet re- 
sistance revealing a large critical regime extending sub- 
stantially above the lowest attained temperature T = 4 
K. From the nonuniversal parameters T c (x) and Rq (x) 
disclosing the respective quantum critical properties we 
derived for the quantum critical exponents the estimates: 
zv ~ 1.46 from T c (x), 77 ~ 0.63 and z ~ 2.35 from 
Ro (x), yielding zv ~ 1.48. Thus, in contrast to the stan- 
dard quantum scaling approach, providing an estimate 
for zv only, the BKT scenario uncovered zv, V and z 
from the quantum critical behavior disclosed in T c (x) 
and Rq (x). Additional evidence for zv ~ 1.5 was estab- 
lished from the comparison of the scaled data with the 
explicit scaling BKT scaling form of the superconductor 
branch. We observed that the scaled data does not fall 
entirely on the BKT curve. Instead a flow to and away 
from the universal characteristics occurred. As T c /T de- 
creases for fixed T c a crossover to AL-behavior sets in, 
while the rounding of the transition leads with increas- 
ing T c /T to a flow away from criticality. The important 
lesson then is that the quality of the data collapse on a 
single curve heavily depends on the temperature range 
of the data entering the plot. Another striking feature is 
the extended scaling regime. Within the BKT scenario it 
follows from the fact that the explicit scaling form of the 
sheet resistance applies along the entire BKT transition 



line irrespective of the distance from the QSI transition. 
In the quantum scaling approach this property remains 
hidden and merely suggests an extended quantum criti- 
cal regime. The piecewise excellent quality of the data 
collapse also reveals that the provision of the substan- 
tial variation of Rq (x) is essential, while in the quantum 
scaling approach it is fixed by the critical sheet resis- 
tance. Supposing that the equivalence between quantum 
phase transitions with D spatial dimensions and the ones 
of classical phase transitions in {D + z) dimensions ap- 
plies, the 2D-QSI transition at the endpoint of a BKT line 
T c (x) should belong to the finite temperature (2 + z)—xy 
universality class. Our estimate z ~ 2.35 and with that 
D = 4.35 — xy critical behavior where 77 = 1/2. However 
the fate of this clean critical point under the influence of 
disorder is controlled by the Harris criterion ! 27 ! 28 If the 
inequality v > 2/D is fulfilled, the disorder does not 
affect the critical behavior. If the Harris criterion is vio- 
lated ( V < 2 /D) , the generic result is a new critical point 
with conventional power law scaling but new exponents 
which fulfill v < 2/D. Since 77 = 1/2 violates this in- 
equality disorder is relevant and drives the system from 
the mean-field to an other critical point with different 
critical exponents as our estimate, 77 ~ 0.63, consistent 
with 77 < 2/D = 1, uncovers. The resulting 2D-QSI 
transition with 77 ~ 0.63 and z ~ 2.35 violates then the 
equivalence between quantum phase transitions in sys- 
tems with D spatial dimensions and the ones of classi- 
cal phase transitions in (D + z) dimensions. In addition 
the proportionality between d/\ 2 (0) and T c (Eq. (HJ)), 
valid below the upper critical dimension D u = 4, does no 
longer hold because (D + z) ~ 4.35 is above D u = 4. In 
fact magnetic penetration depth measurements taken on 
underdoped high quality YBa2Cu3 06+a; single crystals 
revealed T c oc (d/X 2 (0))°' 61 3 Another option is that the 
disorder destroys the QSI transition. Given the evidence 
for a rounded BKT transition (see Fig. [I} and noting 
that the analyzed data do not extend very close to the 
QSI transition (see Fig. [4} further studies are required to 
elucidate this option. In any case, our estimate 77 ~ 0.63, 
the evidence for BKT behavior at finite temperature and 
the Harris theorem imply the presence and relevance of 
disorder at the QSI transition and its irrelevance at fi- 
nite temperature. Further evidence for the importance 
of disorder follows from the temperature dependence of 
the sheet resistance in the insulating phase. Considering 
R (Vg = 0, T) of Bollinger et al— we observe consistency 
with the Mott variable range hopping model in 2D. The 
conductivity exhibits the characteristic temperature de- 
pendence a = 0o exp (^— (T /T) 1 ^ D+1 ^ which applies to 
strongly disordered systems with localized states4^ 

A previous and analogous analysis of the sheet resis- 
tance data of the superconducting LaAlOs/ SrTiOs in- 
terface revealed the critical sheet resistance R c = 2.7 kfi 
and the exponents z77~ 2/3, 77 ~ 2/3 and z ~ 1j2£ In this 
case the gate voltage tuned QSI transition of the 2D elec- 
tron system at the LaAlOa/ SrTiC>3 interface has a finite 
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temperature counterpart in (D + z) dimensions namely 
the (2 + z) = 3-xy model where hyperscaling and with 
that the proportionality between c?/A 2 (0) and T c (Eq. 
(fTTj) ) applies. On the other hand V ~ 2/3 < 2/D = 1 
implies according to Harris theore m 27 i 28 that disorder is 
relevant as well. 

To comment on the BKT-features in the temperature 
dependence of the magnetic penetration depth we con- 
sidered the data of Bert et al££ for a gate voltage tuned 
superconducting LaAlOs/SrTiOa interface in terms of T c 
vs. d/X 2 (T = 0.04 K). The data reveals the flow to the 
universal relationship but much lower T c must be at- 



tained to reach quantum critical regime. However in this 
lowT c regime is the Nelson-Kosterlitz jump given by Eq. 
([TU| very small and appears to be beyond present exper- 
imental resolution^ 2 - Otherwise the data points resemble 
the outline of a fly's wingjil remarkably similar to the 
T c vs. 1/A 2 b (T = 0) plots of the bulk superconductors 
Yo.sCuo.2-123, Tl-1212^ 2 . and Tl-220li£, covering nearly 
the entire doping regime in the so called superconducting 
dome extending from the undcrdoped to the overdoped 
limit. 
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